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ABSTRACT 

Measurements  of  high  frequency  temperature  and  velocity 
fluctuations  were  made  in  the  atmospheric  boundary  layer. 
These  data  were  analyzed  using  analog  correlating  techniques 
to  obtain  the  structure  function  rather  than  spectral  methods 
to  obtain  the  spectral  densities.   The  autocorrelation  func- 
tion of  temperature  fluctuations  at  4  meters  and  2  meters  was 
computed  as  was  the  autocorrelation  function  of  downstream 
velocity  fluctuations  at  4  meters.   The  cross-correlation 
function  of  temperature  and  downstream  velocity  fluctuations 
at  4  meters  was  also  examined.   These  results  supported  the 
Kolmogorov  ideas  at  high  wave  numbers.   There  is  apparently 
very  little  correlation  between  downstream  velocity  fluctua- 
tions and  temperature  fluctuations.   This  indicates  the  inde- 
pendence of  the  two  sensing  systems  and  that  temperature 
fluctuations  are  related  primarily  to  vertical  air  movement 
rather  than  random  fluctuations  in  the  horizontal  air  flow. 
This  aspect  is  further  clarified  by  examining  the  temperature 
fluctuation  records  at  two  heights. 


TABLE  OF  CONTENTS 

I.   INTRODUCTION . 6 

II.   DESCRIPTION  OF  DATA 9 

A.  FIELD  SITE 9 

B.  EXPERIMENTAL  EQUIPMENT  AND  SENSORS  11 

1.  Platinum  Resistance  Thermometer  11 

2.  Hot-Wire  Anemometer  11 

3.  Cup  Anemometer  and  Quartz  Thermometers  12 

4.  Recording  System  12 

III.   THEORETICAL  RELATIONS ■ —  17 

A.  AUTOCORRELATION  AND  CROSS  CORRELATION 
FUNCTIONS 17 

B.  KOLMOGOROV'S  POWER  LAWS  AND  THE  STRUCTURE 
FUNCTION 19 

IV.   ANALYSIS  PROCEDURES  24 

A.  SELECTION  OF  DATA 24 

B.  CORRELATION  EQUIPMENT  AND  PROCEDURES  2  8 

V.   RESULTS 30 

A.  EVIDENCE  OF  ORGANIZED  MOTION  IN  THE  FLOW  30 

B.  CORRELATION  ANALYSIS  34 

C.  STRUCTURE  FUNCTION  40 

VI.   CONCLUSIONS  47 

APPENDIX  A:   Digital  Spectral  Analysis  Procedures  49 

BIBLIOGRAPHY  53 

INITIAL  DISTRIBUTION  LIST  55 

FORM  DD  1473 56 


LIST  OF  FIGURES 

Figure  .       Page 

1.  Location  of  Field  Site 10 

2.  Block  Diagram  of  Recording  System  14 

3.  Data  Record  Selected  for  Analysis  (Direct  and 
Differentiated  Signals)  26 

4.  Data  Record  Selected  for  Analysis  (Direct  Signals)-  27 

5.  Strip  Chart  Record  of  Temperature  Pulses 
Exhibiting  Plume  Properties ' —  32 

6 .  Schematic  of  Thermal  Plumes  in  Shear  Flow  33 

7.  Autocorrelation  of  Direct  Temperature  Signal  at 

4  Meters  for  Various  Lag  Times  (Tau)  35 

8.  Autocorrelation  of  Differentiated  Temperature 

Signal  at  4  Meters  for  Various  Lag  Times  (Tau)  36 

9.  Strip  Chart  and  Autocorrelation  for  Various  Lag 
Times  (Tau)  of  Direct  Temperature  Signal  at 

4  Meters 38 

10.  Cross  Correlation  of  Direct  Temperature  and 
Velocity  Signals  at  4  Meters  for  Various  Lag  Times 

(Tau)  39 

11.  Structure  Function  Plot  of  Direct  Temperature 

Signal  at  2  Meters i*1 

12.  Structure  Function  Plot  of  Direct  Temperature 

Signal  at  4  Meters  42 

13.  Structure  Function  Plot  of  Direct  Velocity  Signal 

at  4  Meters  43 

14.  Structure  Function  Plot  of  Differentiated 
Temperature  Signal  at  2  Meters  

15.  Structure  Function  Plot  of  Differentiated 
Temperature  Signal  at  4  Meters  45 

16.  Structure  Function  Plot  of  Differentiated 

Velocity  Signal  at  4  Meters  *" 


ACKNOWLEDGEMENT 

The  author  gratefully  wishes  to  thank  his  advisor, 
Professor  Noel  E.  J.  Boston,  for  providing  unlimited  hours 
of  his  time  in  helping  the  author  over  the  many  hurdles  and 
hours  of  crises  involved  in  this  thesis.   Without  his  endur- 
ing patience  and  essential  advice,  this  work  would  not  have 
reached  fruition.   The  author  is  also  grateful  for  the  timely 
and  enthusiastic  assistance  of  Professor  Kenneth  L.  Davidson 
whose  encouragement  and  constructive  ideas  were  most  helpful. 

The  author's  thanks  also  go  to  Professor  Thomas  M. 
Houlihan  who  provided  the  electronics  necessary  for  analog 
correlation. 

The  author  would  like  to  acknowledge  the  vital  technical 
assistance  in  computer  programming  provided  by  Mr.  Robert  L. 
Limes  and  Mr.  David  F.  Norman  of  the  Naval  Postgraduate 
School  staff. 

Completion  of  this  project  was  greatly  assisted  by  Lt. 
Robert  D.  Jones  who  generously  shared  his  time  and  knowledge 
in  many  areas  of  data  processing. 


I.   INTRODUCTION 

The  need  to  develop  methods  of  analyzing  high  frequency 
turbulence  signals  is  apparent.   More  emphasis  is  being  given 
to  the  collection  of  high  frequency  turbulence  data  both  in 
the  atmosphere  and  in  the  ocean.   This  is  partly  due  to  the 
fact  that  technological  developments  have  made  such  measure- 
ments feasible  and,  equally  important,  due  to  an  increasing 
need  for  such  data  in  directly  estimating  fundamental  turbu- 
lent parameters  such  as  energy  and  scalar  dissipation  rates. 

Analysis  of  high  frequency  turbulence  data  presents  some 
difficulties  not  encountered  in  low  frequency  analysis.   Ex- 
tremely large  amounts  of  information  are  contained  in  rela- 
tively short  records.   In  many  cases  this  means  the  data  must 
be  analyzed  in  real  time  since  the  compression  of  the  record 
by  playback  speed-up  is  not  usually  possible.   If  analog  data 
are  digitized,  sampling  rates  must  be  very  high  which  limits 
the  length  of  data  record  that  may  be  transferred  to  digital 
tape.   High  frequency  turbulence  signals  tend  to  be  intermit- 
tent so  that  areas  of  strong  correlation,  zero  correlation  or 
low  correlation  may  be  found  over  relatively  short  intervals 
of  two  signals.   Furthermore,  high  frequency  signals  are  small, 
often  not  much  greater  than  system  noise  level. 

These  are  some  of  the  problems  investigated  in  this  the- 
sis.  First,  how  does  one  go  about  analyzing  such  data;  and 
second,  once  analyzed  how  are  results  to  be  interpreted? 


Correlation  and  spectral  analyses  have  found  wide  appli- 
cation in  engineering  and  physical  investigations  for  some 
time.   More  recently  the  techniques  and  experiences  of  these 
established  sciences  have  been  carried  over  into  investigating 
geophysical  phenomena.   Perhaps  the  outstanding  example  in 
meteorology  and  oceanography  of  the  application  of  these  tech- 
niques is  in  micrometeorology  and  air-sea  interaction.   The 
correlation  between  temperature  fluctuations  and  vertical 
velocity  fluctuations  leads  to  a  direct  estimate  of  sensible 
heat  flux,  one  of  the  major  terms  in  thermally-driven  models 
of  the  atmosphere.   The  correlation  between  downstream  and 
vertical  velocity  fluctuations  leads  to  a  direct  estimate  of 
the  Reynolds  stress,  or  momentum  flux,  which  plays  a  dominant 
role  not  only  in  wind-generated  wave  theory  but  also  in  large- 
scale  wind-driven  ocean  circulation.   Spectral  analyses,  com- 
bined with  turbulent  theories,  have  allowed  direct  estimates 
of  the  rate  of  turbulent  kinetic  energy  dissipation  and  the 
rate  of  scalar  dissipation  of  such  quantities  as  temperature 
and  humidity  fluctuations.   Such  estimates  require  measure- 
ments of  turbulent  fluctuations  to  the  highest  possible  res- 
olution of  sensors  and  electronics.   The  distribution  of  these 
quantities  in  space  and  time  leads  to  an  increase  of  under- 
standing of  turbulent  processes.   Whereas  clearly  the  goal  is 
the  understanding  of  the  large  scale  dynamics  of  the  atmos- 
phere and  ocean,  such  understanding  cannot  be  achieved  until 
the  very  small  processes  and  the  non-linear  ways  in  which  they 
enter  into  large  processes  are  understood. 


One  application  of  studies  in  this  field  which  is  of 
special  interest  to  the  military  is  the  possible  improvement 
of  radio  communications  based  upon  an  understanding  of  atmos- 
pheric turbulence.   Electromagnetic  waves  are  refracted  and 
scattered  in  the  atmosphere.   The  nature  of  this  scattering 
and  refraction  is  determined  in  large  part  by  the  local  re- 
fractive index  of  air,  which  in  turn  is  dependent  upon  both 
the  large  and  small  scale  temperature  fluctuations  present. 
Changes  in  the  refractive  index  can  change  the  path  of  radio 
waves  and  seriously  disrupt  both  long  and  short  range  radio 
communications.   Therefore,  increased  understanding  of  tem- 
perature fluctuations  and  the  statistics  of  temperature  fluc- 
tuation distribution,  even  at  very  small  scales,  can  lead  to 
an  improvement  of  communications  by  providing  the  necessary 
information  to  allow  correct  action  in  changing  location  of 
transmitter  or  receiver,  changing  frequency,  or  changing  the 
method  of  propagation  to  suit  atmospheric  conditions. 


II.   DESCRIPTION  OF  DATA 

Based  upon  the  desire  for  information  concerning  the  dis- 
tribution of  temperature  within  a  turbulent  velocity  field 
and  the  possible  correlation  in  such  a  field  between  temper- 
ature and  velocity  fluctuations,  it  was  decided  to  measure 
simultaneous  temperature  and  velocity  signals  with  instrumen- 
tation as  sensitive  as  possible.   The  difficulty  of  accurate- 
ly measuring  temperature  fluctuations  in  a  velocity  field  led 
to  the  development  of  a  platinum  wire  resistance  thermometer 
by  Boston  (1970).   This  instrument  was  used  to  measure  tem- 
perature fluctuations  and  a  hot-wire  anemometer  was  used  to 
measure  velocity  fluctuations.   Measurements  of  temperature 
were  taken  at  4  and  2  meters  above  the  ground,  and  velocity 
was  recorded  at  4  meters  above  the  ground. 

A.   FIELD  SITE 

The  data  examined  in  this  study  were  collected  at  Boundary 
Bay,  British  Columbia  (Figure  1) .   This  location  was  the  third 
of  three  sites  considered.   The  first,  Spanish  Banks,  British 
Columbia  was  not  entirely  satisfactory  because  temperature 
signals  tended  to  be  small  in  the  stable  air  over  the  ocean. 
The  second,  Hanford,  Washington,  was  not  satisfactory  because 
of  the  varied  meteorological  conditions  encountered  in  air 
over  a  desert.   The  third,  Boundary  Bay,  was  a  compromise  be- 
tween the  first  two.   The  air  was  slightly  unstable  as  it 
blew  across  the  beach  off  the  ocean.   The  site  is  a  tidal 
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mud-flat,  which  is  completely  covered  at  high  tide  and  is 
bare  for  about  2 . 5  km  seaward  at  low  tide.   The  site  has  been 
described  in  detail  by  Dobson  (19  69) . 

B.   EXPERIMENTAL  EQUIPMENT  AND  SENSORS 

All  instruments  were  supported  from  a  portable  mast  at 
heights  of  either  4  or  2  meters.   The  signals  from  the  sensors 
were  carried  by  cables  to  a  panel  truck  which  served  as  the 
housing  for  the  recording  equipment.   The  truck  was  30  meters 
downwind  from  the  mast. 

1.  Platinum  Resistance  Thermometer 

The  sensor  of  temperature  fluctuations  was  a  platinum 
resistance  thermometer  consisting  of  a  wire  0.25y  in  diameter 
and  0.25  mm  in  length.   The  wire  formed  one  leg  of  a  Wheat- 
stone  bridge.   A  square  wave  of  constant  amplitude  is  applied 
to  the  bridge  via  an  80  kHz  multivibrator  and  switching  tran- 
sistor.  A  synchronous  detector  is  used  on  the  bridge  output. 
D.C.  amplification  is  performed  by  an  operational  amplifier 
with  fixed  incremental  voltage  gains.   Its  frequency  response 
is  flat  from  D.C.  to  approximately  10  kHz.   The  platinum  re- 
sistance thermometer  electronics  were  placed  at  the  foot  of 
the  mast. 

2.  Hot-wire  Anemometer 

A  constant  temperature  hot-wire  anemometer  system 
(DISA  model  55D05)  was  used.   The  unit  was  battery  operated 
and  used  a  1:1  bridge  ratio  made  with  a  compensating  cable. 
The  hot-wire  probe  consisted  of  a  platinum  wire  5y  in  diam- 
eter and  1.0  mm  in  length.   The  hot-wire  was  mounted 
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approximately  5  cm  above  the  temperature  sensor.   Elec- 
tronically this  system  can  handle  signals  in  the  frequency 
range  from  D.C.  to  100  kHz.   Both  the  time  response  and  wave 
number  limitation  of  the  wire  allow  signals  up  to  1  kHz  to 
be  measured  without  attenuation  at  wind  speeds  of  4  to 
5  m/sec. 

3.  Cup  Anemometer  and  Quartz  Thermometers 

In  addition  to  the  above  instruments,  a  cup  anemom- 
eter and  two  quartz  thermometers  were  used  in  the  experimen- 
tal setup.   The  cup  anemometer,  mounted  on  the  mast  at  4 
meters,  was  used  to  measure  mean  wind  speed  in  order  to  pro- 
vide (1)  an  in-field  check  on  the  hot-wire  anemometer,  and 
(2)  a  mean  wind  speed,  U  (cm/sec) ,  to  apply  in  Taylor's 
hypothesis  in  the  form 

k  =  2_lJL  (1) 

where  k  is  the  radian  wave  number  (cm   )  and  f  is  frequency 
(Hz)  . 

The  quartz  thermometers  were  used  to  get  an  indica- 
tion of  the  relative  stability  of  the  atmosphere.   This  data 
would  be  important  over  a  24-hour  period,  but  not  over  a 
period  as  short  as  8  minutes,  which  was  the  length  of  the 
recorded  data  used  in  this  thesis. 

4 .  Recording  System 

The  signals  from  the  platinum  resistance  thermometer 
and  hot-wire  and  cup  anemometers  were  recorded  on  separate  FM 
(frequency  modulated)  channels  of  a  magnetic  tape  recorder 
(Sangamo  model  3562)  run  at  60  inches  per  second.   The  high 
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frequency  cut-off  of  the  tape  recorder  was  20  kHz.   The 
quartz  thermometer  digital  meter  reading  was  read  onto  the 
voice  edge  track  of  the  tape.   The  recording  system  is  out- 
lined in  Figure  2 . 

Since  both  the  noise  of  the  system  and  the  maximum 
frequency  of  interest  of  the  temperature  signal  were  uncer- 
tain at  the  time  of  the  experiment,  recordings  were  made  both 
without  the  low  pass  filter  (Krohn-Hite  model  3340)  and  with 
it  at  various  settings  to  attenuate  very  high  frequency  noise. 
Data  chosen  to  be  analyzed  in  this  thesis  were  recorded  with  a 
low  pass  filter  cut-off  at  10  kHz  (-3  db  down  frequency) . 

The  signals  from  the  platinum  resistance  thermometer 
(PRT)  and  the  hot-wire  anemometer  were  treated  similarly. 
Each  was  recorded  in  two  different  ways  chosen  to  optimize 
the  incompatible  requirements  of  noise  performance  and  fre- 
quency response.   One  channel  was  unfiltered  (L.F.  GAIN)  and 
the  other  was  differentiated.   The  purpose  of  the  differenti- 
ation circuit  was  to  improve  the  signal-to-noise  ratio  at  high 
frequencies  so  that  the  very  small  scale  fluctuations  could 
be  more  readily  examined.   L.F.  GAIN  provided  gain  not  only 
for  the  direct  signal  but  also  for  the  signal  prior  to  dif- 
ferentiation.  The  high  frequency  gain  control  (H.F.  GAIN) 
determined  the  frequency  at  which  there  was  unity  gain  in  the 
differentiating  circuit.   H.F.  GAIN  was  adjusted  to  provide 
optimum  gain  (prewhitening)  of  the  frequencies  of  interest. 

The  differentiating  circuits  were  matched  linear 
phase  shift  devices  made  from  Philbrick  P85AU  operational 
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amplifiers.   The  circuit  had  a  6  db/octave  gain  up  to  10  kHz 
(depending  on  setting)  and  then  fell  at  6  db/octave. 

Recording  at  60  in/s  with  FM  electronics  ensured  the 
best  signal-to-noise  ratio  and  provided  sufficient' frequency 
response  (up  to  20  kHz)  to  record  faithfully  the  sharp  rises 
and  falls  of  the  differentiated  signals.   It  has  the  disad- 
vantage of  limiting  the  maximum  length  of  a  given  record  to 
8  to  12  minutes,  depending  on  the  length  of  the  tape. 

Specifically  in  the  case  of.  the  data  accumulated  for 
analysis  in  this  thesis,  measurements  of  temperature  and 
velocity  fluctuations  took  place  at  Boundary  Bay  on  the  after- 
noon of  August  8,  1969.   Measurements  were  recorded  on  mag- 
netic tape  No.  20  3(1) .   Record  length  was  8  minutes  commencing 
at  142  7  PDT .   A  platinum  resistance  thermometer  was  located 
on  the  mast  4  meters  above  the  ground.   The  temperature  signal 
was  passed  through  the  PRT  electronics  where  the  signal  was 
amplified,  then  through  the  Krohn-Hite  low  pass  maximum  flat 
filter  with  the  cut-off  frequency  set  at  10  kHz.   The  output 
of  the  Krohn-Hite  filter  was  fed  to  the  Philbrick  amplifier 
manifold  where  both  the  low  frequency  and  differentiated  sig- 
nals were  amplified  again.   The  low  frequency,  undifferen- 
tiated signal  was  then  recorded  on  tape  channel  5.   The  high 
frequency  differentiated  signal  was  recorded  on  channel  7. 

The  DISA  hot-wire  anemometer  also  was  located  on  the 
mast  4  meters  above  the  ground.   The  velocity  fluctuation 
signal  from  this  sensor  was  fed  through  a  Philbrick  amplifier 
manifold  for  signal  gain.   The  low  frequency,  undifferentiated 


15 


signal  was  then  recorded  on  tape  channel  2.   The  high  fre- 
quency, differentiated  signal  was  recorded  on  channel  4. 

In  addition,  at  2  meters  elevation  there  was  another 
platinum  resistance  thermometer.   The  temperature  fluctuation 
signal  from  this  sensor  was  fed  through  and  amplified  by 
separate  PRT  electronics,  and  then  passed  to  a  Krohn-Hite  low 
pass  maximum  flat  filter  with  a  cut-off  frequency  setting  of 
10  kHz.   The  output  of  the  Krohn-Hite  filter  was  then  again 
amplified  by  the  Philbrick  amplifiers.   The  low  frequency, 
undifferentiated  signal  was  then  recorded  on  channel  1,  and 
the  high  frequency,  differentiated  signal  on  channel  3. 

In  the  case  of  all  three  velocity  and  temperature 
measurements,  each  signal  was  amplified  by  the  PRT  electronics 
and  the  Philbrick  operational  amplifiers  to  such  a  degree  that 
complete  advantage  was  taken  of  the  full  dynamic  range  of  the 
tape  recorder. 
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III.   THEORETICAL  RELATIONS 

Correlation  analysis  is  a  method  of  measuring  the  inter- 
dependence of  variables  associated  with  various  processes, 
often  of  a  random  or  periodic  nature.   Simultaneous  time 
functions,  such  as  turbulence  data  measured  at  a  fixed  point 
in  space,  can  be  described  by  correlation  functions.   These 
functions  may  be  of  the  dependence  of  a  variable  at  one  time 
with  itself  at  another  time  (called  the  autocorrelation  func- 
tion) or  the  dependence  of  a  variable  with  another  variable 
(called  the  cross  correlation  function) . 

A.   AUTOCORRELATION  AND  CROSS  CORRELATION  FUNCTIONS 

The  autocorrelation  function  of  a  time  series  g. (t)  at 

the  time  t  and  (t  +  t)  will  be  indicated  by  R. . (t) .   It  is 

obtained  by  averaging  the  products  of  two  values  separated  by 

a  lag  time,  t,  over  an  observation  time,  T.   The  resulting 

average  product  is  a  function  of  x:  strictly 

T 

/  dt 


'ii(T)  =  t^t-  f  9LMg±it+x)dt 


=  gi(t)gi(t+T)  • 


(3) 


At  t  =  0 .  R. . (x)  has  its  maximum  value  which  is  related  to 
the  variance  of  the  signal  by 


Rii<°>  =gi(t)2  -Rii(T»  (4) 


17 


for  all  t .   Clearly 

Rii(T)  =  Rii("T)  •  (5> 

The  form  of  R..(t)  will  be  characteristic  of  the  original 
signal  g. (t) . 

The  autocorrelation  function  is  related  to  the  energy 
spectrum  E..(f)  (sometimes  called  the  power  spectrum)  as 
follows: 


00 

i.  .   (T)  =   |  E.  . 
11  J    11 


R.  .(T)     =  E..(f)     •     COS    2TTfT     •     df  (6) 

o 
and   the   inverse   is: 


oo 


Eii(f)     =    4  R-:-:(T>      '     cos    27TfT     "    dTf  <7) 


which  merely  states  that  R. . (t)  and  E. . (f)  are  Fourier  cosine 
transforms  of  each  other.  One  can  be  found  from  the  other  by 
integration. 

The  theoretical  relations  involved  in  performing  cross 
correlation  of  a  time  series  are  the  same  as  those  for  the 
autocorrelation  function  described  here  with  the  exception 
that  two  different  signals  are  being  correlated  rather  than 
the  same  signal.   That  is,  instead  of  correlating  g. (t)  with 
itself,  g. (t)  is  correlated  with  g.(t). 
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B.   KOLMOGOROV 'S  POWER  LAWS  AND  THE  STRUCTURE  FUNCTION 

Kolmogorov  (19  41)  postulated  that  high  frequency  velocity 
fluctuations  should  be  isotropic  even  though  they  are  gener- 
ated from  and  embedded  in  anisotropic  eddies  of  low  frequen- 
cies.  If  there  is  a  significant  region  of  the  spectrum  at 
high  frequencies  in  which  there  is  no  energy  feeding,  then 
since  there  is  an  energy  cascade  from  lower  frequencies  by 
inertial  transfer,  with  the  pressure  forces  tending  to  trans- 
fer energy  among  components,  it  would  be  reasonable  to  suppose 
that  these  smaller  scales  would  not  exhibit  the  anisotropy 
of  the  low  wave  numbers  and  would  tend  toward  isotropy.   This 
is  known  as  local  isotropy,  and  has  been  supported  by  experi- 
mental evidence  of  Batchelor  (1953) . 

If  the  Reynolds  numbers  are  high  enough,  turbulence  will 
adjust  through  inertial  transfer  and  viscous  dissipation 
until  a  statistical  steady-state,  independent  of  the  initial 
conditions,  is  obtained.   This  steady-state  consists  of  energy 
transfer  into  the  high  frequency  range  that  is  equal  to  energy 
lost  through  viscous  dissipation. 

Kolmogorov  hypothesized  that  at  high  Reynolds  numbers  in 
any  turbulence,  the  small-scale,  high  frequency  eddies  are 
isotropic  and  in  universal  statistical  steady-state.   Their 
properties  should  depend  only  on  the  rate  of  dissipation  of 
energy  (e)  and  on  the  kinematic  viscosity  (v) ,  and  not  direct- 
ly on  time.   Therefore  motion  associated  with  the  steady- 
state  range  is  uniquely  determined  by  the  parameters  e  and 
v  .   Based  on  dimensional  analysis,  Kolmogorov' s  conclusion 
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about  the  form  of  the  three-dimensional  energy  spectrum  was 

=  F(r-)  (8) 

<ev5)1/4 

where  F  is  a  universal  function,  k  is  the  wave  number,  k   is 

s 

the  Kolmogorov  wave  number,  and  1/k   is  the  Kolmogorov 
microscale . 

This  has  been  verified  quite  conclusively  by  Grant, 
Stewart,  and  Moilliet  (1962)  and  Pond,  Stewart,  and  Burling 
(19  6  3) .   From  their  completely  separate  measurements  in  a 
tidal  channel  and  in  the  atmosphere,  all  results  agreed  with 
the  above  functional  relation  over  4.5  decades  of  wave 
numbers . 

Kolmogorov* s  second  hypothesis  concerned  the  inertial  sub- 
range.  If  the  Reynolds  number  is  large  enough,  say,  in  the 
order  of  10  ,  the  energy-containing  eddies  and  the  dissipa- 
tion wave  number  regions  should  be  well  separated  in  wave 
number  space.   There  should,  therefore,  exist  a  range  of  wave 
numbers  on  the  low  wave  number  side  of  the  dissipation  range 
but  at  higher  wave  numbers  than  the  energy-containing  wave 
numbers  in  which  universal  statistical  steady-state  holds, 
but  at  which  there  is  negligible  dissipation  of  energy.   The 
dominant  feature  would  be  an  energy  cascade  by  inertial  forces 
through  this  inertial  subrange.   Since  viscosity  is  not  im- 
portant in  this  wave  number  range,  statistical  properties  are 
determined  by  the  rate  of  dissipation  only.   Thus,  in  the 
inertial  subrange  by  Kolmogorov's  theory,  E.. (k)  should  have 
the  following  functional  relation: 
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E.  .  (k)  =  Ke2/3k"5/3 


where  K  is  the  Kolraogorov  constant.   This  equation  is  empir- 
ical and  is  obtained  by  dimensional  analysis. 

The  results  presented  in  this  thesis  represent  structure 
functions  of  the  temperature  and  velocity  fluctuations  in  the 
atmosphere.   The  structure  function,  D(r),  depends  on  veloc- 
ity or  temperature  differences  between  two  points  in  the 
turbulent  flow,  i.e.,  for  velocity: 


D-  (r)  =   r  /  -l  *     i    \  i2  (10) 

u.        [u(x+r)  -  u(x)J 

where  x  is  a  point  in  the  turbulent  field  and  r  is  the 
separation  distance. 

If  the  turbulence  is  homogeneous,  the  structure  function 
is  related  to  the  correlation  function,  R. . (r) ,  by 

Di±(r)  =  2[RijL(o)  -  Rii(r)].  (11) 

The  inverse  transform  relation  between  the  correlation 
function  and  the  energy  spectrum,  E.^(k),  is 


00 

ill 


Rii(r)  =      Er  (k)  '  cos  kr  *  dk*  (12) 

The  relation  between  the  structure  function  and  the  spectrum 
is  therefore 


00 

III 


Dii(r)  =  2  I  !    (1  "  cos  kr)  #  Eii(k)  '  dk  *       (13) 


For  conditions  of  isotropy,  equation  (13)  reduces  to 
(Strohbehn  19  6  8) : 
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00 

,  f  , ,    sin 


/  (1  _  —  kr)  .  E..(k)  •  k2  •  dk  .       (14) 
H        I        kr       n 

o 

It  can  be  seen  from  equation  (14)  that  the  structure  func- 
tion, as  the  spectrum  E. . (k) ,  should  also  take  on  a  charac- 
teristic power  law  behavior  in  isotropic  turbulence.   From 
dimensional  reasoning  and  the  assumption  of  dependence  on 
£  and  r,  this  power  law  relation  is 

Du(r)  =  C  e2/3  r2/3  ,  (15) 

where  C  is  a  constant  of  proportionality. 

Similar  arguments  may  be  applied  to  temperature  fluc- 
tuations to  yield  expressions  for  the  temperature  spectrum  and 
temperature  structure  function  with  the  same  power  law  be- 
havior.  For  isotropic  turbulence,  the  temperature  structure 
function  can  be  related  to  the  temperature  autocorrelation 
function  by 


(T1  -  T2)2  -  2T2(1  -  R1±)  (16) 

where  R. .  is  the  autocorrelation  function.   Thus 
n 

(1  -  Ri±)  a  r2/3  .  (17) 

This  relation  is  expressed  in  terms  of  the  separation 
distance,  r,  whereas  the  correlation  analysis  was  made  by 
varying  lag  time,  x.   Assuming  homogeneity  of  the  turbulence 
field,  which  seems  entirely  justified  for  the  observations 
made  in  this  study,  the  correlation  in  space  should  be  equiv- 
alent to  the  correlation  in  time.   The  relation  between  r  and 
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T  IS 

T  =  §■  (18) 

where  x  is  the  lag  time  and  U  is  the  mean  wind  speed.   This 
is  a  variation  of  Taylor's  frozen  turbulence  hypothesis, 
which  .is  most  often  used  to  transform  frequency  (f)  (time) 
spectra  to  wave  number  (k)  (space)  spectra: 

k  =  !-§-£  .  d9) 

From  (18)  and  (19) ,    one  may  relate  t  to  frequency  since 

-  -  k  .  (20) 

r 

So   f  =  ^ .  (21) 

2  If  T 
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IV.   ANALYSIS  PROCEDURES 

Because  of  the  need  to  develop  methods  of  analyzing  high 
frequency  turbulence  signals,  it  was  decided  to  use  two  dis- 
tinct analysis  procedures.   An  attempt  was  first  made  to 
perform  spectral  analysis  on  a  6 -minute  record  of  data  by 
digital  procedures.   This  was  accomplished  by  means  of  analog- 
to-digital  conversion  (XDS  9  300/CI  5000  Hybrid  Computer  System) 
followed  by  IBM-360/6  7  digital  computer  spectral  analysis. 
This  procedure  met  with  only  limited  success  and  is  explained 
in  Appendix  A.   The  second  procedure  consisted  of  a  correla- 
tion analysis  of  the  analog  temperature  and  velocity  fluctua- 
tion signals  by  a  DISA  Analog  Correlator.   The  results  of  this 
second  method  of  analysis  comprise  the  major  results  of  this 
thesis. 

A.   SELECTION  OF  DATA 

Since  the  objectives  of  this  thesis  were  concerned  with 
the  examination  of  high  frequency  turbulence  signals,  it  was 
not  necessary  to  examine  long  sections  of  data.   Thus  it  was 
decided  that  a  record  of  about  one  minute  length  would  be 
sufficient.   Actually  the  exact  record  length  turned  out  to 
be  50  seconds  because  of  the  scale  used  in  recording  the  out- 
put of  the  DISA  Analog  Correlator. 

The  primary  criterion  for  selecting  a  section  of  data  for 
analysis  was  stationarity  of  velocity  signals.   In  addition 
it  was  desirable  that  the  selected  section  contain 
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temperature  fluctuations  at  both  4  and  2  meters  elevation  that 
could  be  considered  typical  of  the  8-minute  record.   There- 
fore, with  these  ideas  in  mind,  the  data  record  selected  was 
a  50-second  section  of  tape  20  3(1)  commencing  22  seconds 
after  the  beginning  of  the  measurements.   This  section  ap- 
pears in  Figure  3.   Note  that  high  frequency  turbulence  inten- 
sity of  both  the  direct  and  differentiated  velocity  signals 
remains  reasonably  uniform  over  the  length  of  this  section. 
Figure  4  shows  the  direct  signals  only.   The  direct  velocity 
signal  shows  that  some  low  frequency  oscillations  are  present, 
but  it  is  not  difficult  to  see  that  the  velocity  signal  devi- 
ates very  little  from  the  mean  wind  speed  of  3.0  m/sec. 

A  variety  of  temperature  fluctuations  were  encountered 
throughout  the  8-minute  record.   Some  sections  showed  a  signal 
that  tended  to  be  one-sided  but  with  the  small  scale  fluctua- 
tions remaining  uniformly  active  throughout.   The  one-sided 
feature  is  a  shape  which  is  often  characteristic  of  tempera- 
ture signals;  a  gradually  rising  leading  edge  and  a  sharp 
drop-off  of  the  trailing  edge,  which  gives  the  record  a  saw- 
tooth pattern.   In  general  the  signal  exhibited  intermittent 
quiet  and  active  periods.   The  50-second  record  selected  at- 
tempted to  incorporate  as  many  of  these  fluctuation  character- 
istics as  possible.   These  characteristics  will  be  interpreted 
in  more  detail  in  view  of  the  scope  and  objectives  of  this 
thesis  in  Section  V. 
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Figure    3.      Data   Record   Selected    for   Analysis 

(Direct   and   Differentiated    Signals) 
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B.   CORRELATION  EQUIPMENT  AND  PROCEDURES 

The  correlation  analysis  was  performed  using  a  DISA  Type 
55D70  Analog  Correlator  System.   In  support  of  the  DISA  cor- 
relator were  a  DISA  Type  55D75  Time  Delay  Unit  and  a  DISA 
Type  52B01  Sweep  Drive  Unit.   The  sweep  drive  unit  was  used 
to  provide  a  fractional  breakdown  of  the  delay  range  setting 
of  the  time  delay  unit.   This  feature  permitted  the  use  of 
convenient  fractional  lag  times  within  the  range  of  the  set- 
ting on  the  time  delay  unit.   In  addition  the  sweep  drive 
unit  provided  control  for  the  Hewlett/Packard  Moseley  135  X-Y 
Recorder  which  was  used  to  record  the  correlator  output. 

The  signals  being  correlated  were  reproduced  on  the  Sangamo 
3562  tape  recorder  and  fed  to  the  time  delay  unit.   The  speed 
of  the  tape  recorder  for  all  correlation  measurements  was 
60  in/s  (real  time) .   The  delay  range  setting  of  the  time 
delay  unit  was  0-10  ms  when  correlating  undifferentiated 
signals  and  0-1  ms  when  autocorrelating  differentiated  signals. 

The  signal  was  then  fed  through  the  DISA  correlator  and 
sweep  drive  unit.   The  integrating  time  constant  setting  used 
in  all  correlations  was  10  seconds.   This  provided  enough 
signal  smoothing  to  recognize  the  large  fluctuations  in  cor- 
relation, while  also  smoothing  the  signal  of  small  scale 
fluctuations  irrelevant  to  the  fixed-lag  correlation  study. 
The  lag  times  set  into  the  sweep  drive  unit  varied  with  the 
correlation  of  different  signals.   The  sweep  drive  unit  had 
the  capability  of  reducing  the  delay  range  setting  into  frac- 
tions of  1/500 .   Thus  different  settings  of  lag  time  could 
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be  chosen  to  adequately  cover  the  correlation  range  from 
-1.0  to  1.0.   Some  specific  settings  for  selected  autocor- 
relations are  given  in  Section  V. 
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V.   RESULTS 

The  results  of  this  study  are  significant  in  light  of  its 
objectives,  which  were  to  (1)  investigate  methods  of  examin- 
ing universal  properties  of  high  frequency  turbulence  signals 
and  (2)  to  specifically  examine  the  relation  between  tempera- 
ture fluctuations  and  downstream  velocity  fluctuations  in  the 
atmospheric  boundary  layer.   A  visual  analysis  of  the  tempera- 
ture fluctuations,  especially  at  4  meters  elevation,  provided 
substantial  support  of  the  work  of  Taylor  (19  5  8)  and  Priestly 
(1959).   The  correlation  analysis  of  both  direct  and  differen- 
tiated signals  at  2  and  4  meters  elevation  obeyed  correlation 
analysis  theory  and  also  agreed  with  Kolmogorov's  -5/3  power 
law  relation  for  the  inertial  subrange  when  the  data  were 
analyzed  by  the  structure  function.   Finally,  a  lack  of  cor- 
relation between  velocity  and  temperature  fluctuations  at  the 
same  height  indicated  the  independence  of  the  two  sensing  sys- 
tems ,  and  therefore  the  fact  that  their  design  is  good. 

A.   EVIDENCE  OF  ORGANIZED  MOTION  IN  THE  FLOW 

A  consistent  characteristic  of  the  temperature  data  was 
the  occurrence  of  "pulses"  of  large  amplitude.   These  pulses 
occurred  nearly  simultaneously  at  both  2  and  4  meters  eleva- 
tion.  Their  occurrence  was  intermittent  and  always  repre- 
sented an  increase  in  temperature.   They  occurred  within 
quiescent  periods  during  which  the  recorded  trace  appeared  to 
have  a  preferred  mean  value. 
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A  selected  portion  of  a  strip  chart  record  exhibiting 
these  pulses  appears  in  Figure  5.   The  pulses  (e.g.,  point  A 
in  Figure  5)  have  a  triangular  shape  which  is  distorted  to 
give  a  saw-tooth  pattern.   Taylor  (195  8) ,  Webb  (19  65)  and 
Priestly  (1959)  have  reported  similar  occurrences  in  data  of 
these  kind. 

Taylor  (195  8)  argues  convincingly  that  the  pulses  are  due 
to  organized  thermal  structures,  "plumes,"  embedded  within 
the  turbulent  flow.   Taylor  reached  his  conclusion  by  perform- 
ing correlation  analyses  between  temperature  fluctuations  at 
several  levels.   The  time  lag  for  maximum  correlations  was  in 
agreement  with  an  expected  tilt  of  these  plumes  due  to  the 
shear  flow.   A  schematic  of  the  plume  in  a  shear  flow  appears 
in  Figure  6 . 

A  time  lag  between  the  occurrence  of  the  pulse  at  the  2 
and  4  meter  level  is  evident  at  Point  A  in  Figure  5.   This 
lag  indicates  that  the  plume  was  swept  past  the  4  meter  level 
before  the  2  meter  level.   Because  the  sensors  were  mounted  in 
a  plum-line  arrangement,  these  results  are  in  agreement  with 
the  concept  that  plumes,  under  the  influence  of  the  wind 
shear,  are  responsible  for  the  observed  pulses. 

The  saw-tooth  characteristic  of  the  pulses  is,  perhaps, 
due  to  the  plumes  having  a  well-defined  trailing  edge,  but  a 
leading  edge  that  is  altered  by  mixing  or  entrainment.   Mixing 
would  lower  the  temperature  in  the  leading  portion  of  the 
plume  giving  the  saw-tooth  pattern. 
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Considering  the  consistent  occurrence  of  pulses  in  these 
temperature  data  in  view  of  Taylor's  (1958)  more  extensive 
study,  it  is  evident  that  organized  motion  exists  within  the 
turbulent  flow.   By  selecting  a  coordinate  system  consisting 
only  of  a  trace  of  the  temperature  versus  time,  it  is  pos- 
sible to  obtain  a  description  of  the  plumes  in  the  complex 
regime.   The  presence  of  organized  motion  does  not  appear  to 
influence  the  universal  properties  of  the  fluctuations  at  the 
small  scales.   Properties  at  high  frequencies  are  examined  by 
another  coordinate  system,  the  structure  function  versus  time 
lag. 

B.   CORRELATION  ANALYSIS 

The  results  of  the  autocorrelation  analysis  were  consist- 
ent and  as  expected.   The  correlation  function,  R. . ,  of  both 
the  direct  and  differentiated  signals  decreased  as  the  lag 
time,  t ,  increased.   This  was  true  throughout  the  50 -second 
record.   Figure  7  is  a  typical  example  of  such  behavior. 
Here  direct  temperature  fluctuations  at  4  meters  elevation 
are  examined.   Note  that  with  x  set  at  0  ms ,  the  autocorrela- 
tion function  exhibits  a  high  value  (-.9  3)  and  that  as  the 
laq  time  is  increased,  R. .  decreases  in  value  until,  with  a 
lag  time  of  10  ms ,  a  negative  correlation  of  mean  value  -.2  8 
exists.   The  fact  that  a  correlation  of  1.0  does  not  result 
from  a  x  setting  of  0  ms  is  attributed  to  correlation  equip- 
ment calibration.   When  the  time  delay  unit  was  set  for  0  ms , 
in  fact  a  slight  time  lag  was  present. 
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Figure  8  exhibits  an  example  of  the  autocorrelation 
results  for  a  differentiated  signal.   In  this  case,  the  dif- 
ferentiated temperature  at  4  meters  elevation  was  correlated 
with  itself.   Since  in  differentiated  signals  the  very  high 
frequencies  (i.e.,  very  small  scale  fluctuations)  are  accen- 
tuated, the  lag  times  used  for  differentiated  signals  are  one 
order  of  magnitude  smaller  than  those  used  in  correlation 
analysis  of  the  non-differentiated,  direct  signals.   However, 
the  same  autocorrelation  characteristics  appear.   When  t  is 
set  at  0 ,  R. .  is  high,  and  as  the  lag  time  is  increased,  the 
value  of  the  correlation  function  decreases. 

In  Figures  7  and  8,  a  definite  pattern  of  autocorrelation 
function  fluctuation  exists  across  the  entire  50-second  record, 
Times  of  interesting  activity  were  selected  and  the  structure 
function  plotted. 

Figure  9  shows  the  autocorrelation  function  of  temperature 
fluctuations  at  4  meters  elevation  matched  with  the  50-second 
analog  trace  of  these  same  temperature  fluctuations.   It  can 
be  seen  that  the  area  of  the  temperature  plume  at  the  36.5 
second  point  coincides  with  a  definite  increase  in  correlation 
at  the  same  time  in  the  autocorrelation  record. 

Figure  10  shows  a  plot  of  the  cross  correlation  of  direct 
signals  of  temperature  and  velocity  fluctuations  at  4  meters 
elevation.   The  continuous  low  values  of  the  correlation 
function,  R.  .,  regardless  of  the  lag  time,  x,  indicates  that 
these  two  signals  are  poorly  correlated.   This  lack  of  cor- 
relation between  simultaneous  temperature  and  velocity 
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fluctuations  at  the  same  height  indicates  that  the  design  of 
the  two  sensors,  the  platinum  resistance  thermometer  and  the 
hot-wire  anemometer,  is  favorable  for  measuring  such  signals, 
A  lack  of  correlation  indicates  that  velocity  fluctuations 
are  not  having  a  significant  effect  on  the  platinum  resis- 
tance thermometer;  and  vice  versa,  that  the  temperature  fluc- 
tuations do  not  affect  the  velocity  measurements  of  the  hot- 
wire anemometer. 

C.   STRUCTURE  FUNCTION 

Figures  11  through  16  show  logarithmic  plots  of  the  auto- 
correlation structure  function,  1-R. . ,  versus  lag  time,  t, 
for  all  six  recorded  signals  at  various  times  of  interest  in 
the  50 -second  record.   In  the  cases  of  direct  temperature 
fluctuations  at  2  and  4  meters  elevation  and  direct  velocity 
fluctuations  at  4  meters  elevation,  the  slope  of  the  plots 
approximates  the  +2/3  value  predicted  by  applying 
Kolmogorov's  inertial  subrange  hypothesis  to  the  structure 
function.   The  slopes  of  the  differentiated  signals  decrease 
at  large  values  of  t . 
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VI.   CONCLUSIONS 


Analysis  of  the  strip  chart  records  of  temperature  fluc- 
tuations indicated  that  the  data  support  Taylor's  and 
Priestly *s  description  of  thermal  plumes  within  a  turbulent 
shear  flow;  and,  as  such,  this  study  added  information  on 
these  phenomena. 

The  data  were  consistent  with  predictions  by  similarity 
theory  since  both  velocity  and  temperature  fluctuations 
assumed  a  +  2/3  power  law  behavior  for  the  structure  function 
in  accordance  with  Kolmogorov's  inertial  subrange  hypotheses. 

The  results  of  this  study  provide  further  support  that 
the  +2/3  power  law  behavior  is  valid  for  temperature  fluc- 
tuations (Krechmer,  1952) .   Results  of  previous  work  in  the 
atmospheric  boundary  layer  have  been  primarily  concerned  with 
the  structure  function  behavior  of  velocity  fluctuations.   At 
first  glance,  it  is  surprising  that  temperature  as  well  as 
velocity  fluctuations  should  exhibit  consistent  +2/3  slopes 
since  the  temperature  traces  showed  the  presence  of  intermit- 
tent active  and  quiet  periods  associated  with  thermal  plumes. 
These  plumes  should  be  associated  with  anisotropic  conditions 
since  they  imply  vertically  upward  directional  preference  due 
to  buoyancy.   Thus,  the  isotropy  deduced  from  the  +2/3  behav- 
ior of  the  structure  function  and  the  anisotropy  of  the 
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characteristic  plumes  appear  to  be  a  paradox.    This  apparent 
inconsistency  can  be  explained  by  the  fact  that  in  this  study 
only  the  very  high  frequency  temperature  fluctuations,  which 
are  independent  of  the  features  of  the  large  scale  flow, 
were  examined.   In  order  to  see  the  difference  in  spectra  or 
structure  functions,  the  lower  frequencies  must  be  examined. 
The  examination  of  these  lower  frequencies  is  one  obvious 
future  extension  of  this  work. 

The  favorable  design  of  the  two  sensors  as  indicated  by 
the  lack  of  correlation  between  velocity  and  temperature 
fluctuations  at  4  meters  elevation  suggests  that  the  two  sens' 
ing  systems  used  to  measure  the  turbulence  in  this  study  are 
compatible  and  could  continue  to  be  used  to  make  further 
measurements . 


A  +2/3  slope  does  not  guarantee  isotropy;  isotropy  is 
a  sufficient  condition  for  a  +2/3  slope,  but  not  a  necessary 
condition.   For  examples,  see  Weiler  and  Burling  (1967). 
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APPENDIX  A 

DIGITAL  SPECTRAL  ANALYSIS  PROCEDURES 

A  significant  portion  of  the  research  in  this  study  was 
an  attempt  to  digitize  the  analog  signals  of  temperature  and 
velocity  fluctuations,  and  then  perform  digital  spectral 
analysis  upon  the  data  by  means  of  a  high  speed  digital  com- 
puter.  Despite  the  fact  that  various  problems  were  encoun- 
tered and  only  limited  success  resulted,  it  is  worthwhile  to 
discuss  the  procedures  involved  in  this  analysis  since  they 
could  well  be  applicable  to  future  studies  of  the  relation 
between  velocity  and  temperature  fluctuations. 

A.   ANALOG  TO  DIGITAL  CONVERSION 

The  analog  to  digital  conversion  was  accomplished  by 
means  of  a  XDS  9  300/CI  5000  Hybrid  Computer  system.   The  XDS 
9  300  is  a  digital  computer  manufactured  by  Xerox,  Inc.  with 
a  32,000  byte  storage  capacity,  and  the  CI  5000  is  a  COMCOR 
analog  computer  with  a  50-amplifier  size  rating.   The  temper- 
ature/velocity signal  output  of  the  Sangamo  3562  tape  recorder 
was  connected  to  a  Krohn-Hite  model  3340  low  pass  filter.   The 
cut-off  frequency  settings  were  2  kHz  for  the  undifferentiated 
signals  and  4  kHz  for  the  differentiated  signals.   The  signals 
were  then  fed  to  the  analog  computer  patch  panel.   The  signal 
was  amplified  in  this  patch  panel  to  take  advantage  of  the 
full  dynamic  range  of  the  hybrid  computer  system.   In 
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subsequent  troubleshooting,  the  patch  panel  was  found  to  be  a 
source  of  substantial  noise.   It  is  possible  that  this  was 
one  of  the  problems  that  precluded  good  results  from  this 
procedure . 

From  this  patch  panel,  the  signals  were  fed  to  the  XDS 
9  300  digital  computer.   The  digitized  output  of  the  XDS  9  300 
computer  was  then  recorded  on  7 -track  magnetic  tape.   Two 
channels  of  the  analog  tape  recorder  were  digitized  at  one 
time.   The  procedures  resulted  in  digital  tape  containing 
blocks  of  data  2048  words  in  size.   Every  other  word  iri  each 
block  was  a  data  sample  from  one  channel.   Thus  each  block  of 
data  contained  10  24  words  (data  samples)  of  each  channel. 
The  sampling  rate  used  for  digitizing  the  undifferentiated 
signal  was  2000  samples/second.   The  differentiated  signal 
was  sampled  at  a  rate  of  4000  samples/second.   In  the  case  of 
all  signals,  a  6-minute  section  of  the  8-minute  record  was 
digitized.   The  procedure  of  digitizing  two  channels  at  one 
time  was  another  possible  source  of  error  in  this  process. 
In  future  work,  it  is  recommended  that  initially  only  one 
channel  be  digitized  at  a  time  in  order  to  simplify  trouble- 
shooting in  later  stages  of  this  procedure.   Of  course,  when 
computing  cross-spectra,  it  will  be  necessary  to  digitize 
two  channels  simultaneously. 

For  a  more  specific,  detailed  description  of  the  general 
analog  to  digital  conversion  procedure  used  in  this  study, 
including  programs,  equipment  settings,  and  further  recom- 
mendations, see  Jones  (19  71) 
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B.   DIGITAL  DATA  PROCESSING 

The  digital  data  processing  was  performed  on  an  IBM 
360/6  7  computer.   This  system  contains  a  core  space  of 
762 , 000  bytes.   Supporting  equipment  includes  two  compilers, 
two  CALCOMP  plotters,  two  printers,  and  twenty  remote  units. 
Digital  tapes  created  by  the  hybrid  computer  system  were 
written  in  7-track  form  (octal  representation) .   However  the 
IBM  360/6  7  System  normally  operates  in  9 -track  form  (hexa- 
decimal representation) .   Therefore  all  data  derived  from  the 
analog  to  digital  conversion  had  to  be  converted  to  a  9 -track 
format  prior  to  analysis  by  the  IBM  360/6  7.   Standard  sub- 
routines existed  for  this  conversion,  but  required  some  modi- 
fications as  explained  in  Jones  (19  71) .   The  output  of  this 
conversion  was  the  digitized  data  on  a  9 -track  tape.   This 
process  was  another  possible  source  of  error.   In  trouble- 
shooting this  phase  of  the  procedure,  a  program  was  written 
to  plot  the  data  of  the  9 -track  tape.   Numerous  spikes  not 
present  in  the  raw  analog  data  were  found  to  be  present  in 
the  9 -track  data. 

The  next  step  in  this  procedure  was  the  computation  of 
the  Fourier  coefficients  for  the  data.   The  program  FTOR, 
described  by  Wilson,  Boston,  and  Denner  (1969),  was  used  to 
accomplish  this  step.   The  output  of  the  Fourier  coefficients 
was  recorded  on  another  9 -track  tape. 

The  final  step  in  the  digital  data  processing  was  the 
computation  and  plotting  of  the  spectra.   Two  plotting  pro- 
grams, SCOR  and  FCPLOT,  described  by  Wilson  et.  al .  (1969), 
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computation  and  plotting  of  the  spectra.   Two  plotting  pro- 
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read  the  output  tape  of  the  FTOR  program  and  plotted  various 
spectra.   The  SCOR  program  provided  a  plotted  output  of 
power  spectra  and  cross  spectra,  and  the  FCPLOT  program  pro- 
duced printer  and  ink.  drawings  of  amplitudes  of  Fourier  coef- 
ficients.  Both  programs  processed  a  single  reel  of  digital 
tape  (1368  blocks  of  2048  words  each)  in  less  than  10  minutes 
computing  time.   In  this  study  the  SCOR  program  output  of 
spectra  and  cross  spectra  for  all  measured  temperature  and 
velocity  fluctuations  exhibited  little  relation  to  known  power 
spectra  behavior  for  such  signals.   The  FTOR,  SCOR,  and  FCPLOT 
programs  have  been  used  successfully  in  the  past  by  Wilson, 
et.  al^.  (1969)  in  analyzing  such  data.   This  leads  to  the  be- 
lief that  the  problems  encountered  and  general  lack  of  results 
in  this  study  were  probably  generated  in  either  the  analog  to 
digital  conversion  or  the  7-track  to  9 -track  tape  conversion 
processes . 

Detailed  information  on  the  FTOR,  SCOR,  and  FCPLOT  pro- 
grams can  be  found  in  Wilson,  et.  al_.  (1969)  .   Specific  de- 
tails of  the  application  of  these  programs  to  temperature  and 
velocity  fluctuations  in  the  atmospheric  boundary  layer  are 
described  at  length  by  Jones  (19  71). 
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